Draft version July 11, 2012 

Preprint typeset using I^'T^]X style emulatcapj v. 12/16/11 



FIRST-2MASS RED QUASARS: TRANSITIONAL OBJECTS EMERGING FROM THE DUST 

EiLAT Glikman^ Tanya Urrutia^, Mark Lacy^, S. George Djorgovski"*, Ashish Mahabal^, Adam D. Myers^, 
Nicholas P. Ross'*, Patrick Petitjean", Jian Ge**, Donald P. Schneider'-' "', Donald G. York " 

Draft version July 11, 2012 

ABSTRACT 

We present a sample of 120 dust-reddened quasars identified by matching radio sources detected at 
1.4 GHz in the FIRST survey with the near-infrared 2MASS catalog and color-selecting red sources. 
Optical and/or near-infrared spectroscopy provide broad wavelength sampling of their spectral energy 
distributions that we use to determine their reddening, characterized hy E{B — V). We demonstrate 
that the reddening in these quasars is best-described by SMC-like dust. This sample spans a wide 
range in redshift and reddening (0.1 z 0.1 < E{B — V)< 1.5), which we use to investigate 
the possible correlation of luminosity with reddening. At every redshift, dust-reddened quasars are 
intrinsically the most luminous quasars. We interpret this result in the context of merger-driven 
quasar/galaxy co-evolution where these reddened quasars are revealing an emergent phase during 
which the heavily obscured quasar is shedding its cocoon of dust prior to becoming a "normal" blue 
quasar. When correcting for extinction, we find that, depending on how the parent population is 
defined, these red quasars make up < 15 — 20% of the luminous quasar population. We estimate, 
based on the fraction of objects in this phase, that its duration is 15 — 20% as long as the unobscured, 
blue quasar phase. 

Subject headings: dust, extinction; quasars: general; surveys 



1. INTRODUCTION 

Quasars are among the most energetically powerful 
objects in the universe. Their luminosities are indica- 
tive of how supermassive black holes (SMBHs; 10® — 10^ 
Mq) grow; namely, through accretion at or near the Ed- 
dington limit. Many observations suggest that SMBH 
growth and the build-up of their host galaxies are closely 
linked. SMBHs appear to be a ubiquitous feature at 
the centers of all large galaxies. SMBH masses are pro- 
portional to the mass/velocity dispersio n of their host 
spheroid (the M— cr relation: .Ferrare sc fc M errittll200fll : 
iGebhardt et all 120001 : iMagorrian et'a l. 1998). Further- 
more, the cosmological evolution of the QSO luminos- 
ity density a nd of the star formation rate are r emark- 
ably similar f Shaver et al.1[T996l: I Wall et al.l[2005h . The- 
oretical models of cosmological galaxy evolution require 
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feedback from the SMBH growth as an ingredient in or- 
der to reproduce the observed fiT-band number coiuits 
and g alaxy luminosity function (iKauffmann &: Haehneltl 
[200g ^ lDi Matteo et al.ll2005HCroton et al.ll2006D . 

All of this information forms the basis for quasar- 
galaxy CO- formation models. In these models, SMBHs 
and galaxies co-evolve through major mergers which fuel 
both a starburst and acc r etion onto the nuclea r black 
hole (jSanders et al.l 11981 iHopkins et all l2005aL [2006). 
The models predict an obscured phase for the young 
quasar resulting from the large amounts of gas and dust 
tunneled inwards to fuel the quasar. A luminous blue 
quasar emerges at the end of this process, shining until 
it exhausts its fuel supply. In this picture reddening is 
correlated with the evolutionary stage of the quasar. 

Such scenarios have been invoked to explain the pres- 
ence of buried AGN seen in ultraluminous infrared 
galaxies (ULIRGSs; Sanders et a l. 1988), a high frac- 
tion of w hich also show evidence of merging and in- 
teraction (jSanders fc Mirabell Il996[ ). However, lumi- 
nous, blue quasars show few signs of inte raction; their 
hosts are mostly undisturbed galaxies fe.g..'D unlop et aTl 
.2003; Flovd et al. 2004; Zakamska et al. 200i, but see 
Bennert et al. 2008). These seemingly conflicting obser- 
vations suggest that dust that obscures buried quasars 
might be cleared during a transitional phase. Such tran- 
sitional objects would represent a key link in the evolu- 
tionary path, where the dust that completely obscured 
the AGN hosted in the ULIRG is cleared, as signs of in- 
teraction dissipate, eventually to reveal an unobscured, 
luminous quasar. This missing link should be a popu- 
lation of highly reddened, but not completely obscured, 
quasars which have been largely missed by optical quasar 
surveys. 

Obscuration in AGN is often discussed in the con- 
text of an orientation-based AGN-unification model (c.f., 



Antonuccil [19931: lUrrv fc Padovanil 119951: iTreister et all 
2004f ). This model provides an inherent description of ac- 



creting black holes as having an axisymmetric geometry 
with an accretion disk surrounded by photoionized clouds 
that give rise to broad emission hues. Dusty, molecular 
clouds, often modeled as a torus co-planar with the ac- 
cret ion disk (thou gh the structure may be more complex, 
e.g. ■ lElitzuill2008f ) , obscure the broad emission lines along 
certain lines-of-sight. Farther out lie cooler photoionized 
clouds emitting narrow forbidden lines. AGN seen along 
certain lines of sight intersect the obscuring torus and 
reveal only narrow-emission lines in their spectra. 

The unification scheme is well-established for low- 
luminosity AGN, and has been extended to so-called 
Type 2 (obs cured) quasars as well (jZakamska et al.|[2003l 
12004 120051 ). While Type 2 quasars are interesting in 
their own right, recent research confirms that these ob- 
jects are not reddened as a result of recent merger and 
are not tied into the evolutionary picture of galaxy and 
quasa r co-formation ([Sturm et al.ll2006l : iZakamska et al.l 
l2006f ). Observations of low- luminosity AGN selected in 
X-rays or the mid-infrared show a strong increase in the 
ratio of obscured/unobscured AGN with decreasing lu- 
minosity (e.g. jTreister et al.ll2008l . l2009[) r'l One possible 
interpretation of this is that the covering fraction from 
the obscuring material shrinks with inc reasing AGN lu - 
minosity (the so-called "receding torus" lLawrenc"3l 1 99 ll ) . 

In this paper, we make a distinction between reddened 
Type 1 objects (whose spectra show broad-line emission) 
and Type 2 AGN which have only narrow emission lines 
and are typically more heavily obscured. We demon- 
strate that this distinction selects a unique population 
of objects whose physical mechanism for the reddening 
appears to differ from an orientation-based extinction. 

Reddened quasars can be missed for several reasons: 
(1) dust extinction can dim sources below the detection 
limit of a given imaging survey, and (2) dust-reddening 
changes the color of a quasar, reddening ultraviolet ex- 
cess (so-called UVX) selected quasars into the optical 
color locus of stars as a complex function of redden- 
ing and redshift - particularly at redshifts of around 
1.5 < z < 2.5, where quasars are most abundant. (3) 
Morphological selection of point sources excludes quasars 
whose rest-frame UV and optical emission is dimmed 
relative to their host galaxies, causing them to appear 
extended, especially in blue-band images. To mitigate 
these obstacles, red quasars must be selected at wave- 
lengths that are less sensitive to dust extinction (e.g., 
radio, infrared, or X-rays), using color cuts that include 
dust-reddened objects. 

Small samples of red quasars have been found using 
various methods - selection in the radio (Webster et al 



19951 iGregg et al.| [2b02: Wh ite et al.1' 2003: Urrutia et al 
2009), near-infrared (Cutri et al.l 12001; Maddox et al 



2008; Gcorgakakis et al. 2009), mid-infrared (jLacv et al 
J004, 2007; PoUctta ct al. j2006l i2008j), and X-rays 
{Kim & Elvis 1999; Pollett a et al.l I2OO6I ) - with large 
variations in the estimated fraction of red quasars in 
the overall quasar population. If red quasars represent 
the transitional phase between hidden accretion and un- 
obscured radiation, then the relative abundance of red 



quasars constrains the duration of this transitional phase. 
Furthermore, catching quasars in the act of this "clearing 
out" may inform the underlying processes linking galaxy 
evolution to black hole growth. 

In this paper we present a complete sample of 120 dust- 
reddened quasars, with well-understood selection criteria 
~ the largest to date. We use this sample to constrain the 
fraction of reddened quasars. We also study their red- 
dening properties and the evolution of these properties 
with redshift and luminosity. 

The paper is organized as follows: we review the results 
from our previous work and motivate our color selection 
for finding red quasars in Section [5J We describe the 
candidate selection process as well as followup spectro- 
scopic observations and classification in Section |3l In 
Section |4l we describe the properties and demograph- 
ics of this red quasar population compared to optically- 
selected quasars and in particular UVX-selected quasars. 
In Section [5l we explore the nature of the reddening and 
its effect on estimating the fraction of reddened quasars. 
Section [6] concerns a unique subset of broad absorption 
line (BAL) quasars whose fraction is significantly higher 
among red quasars than in the overall quasar popula- 
tion. Finally, in Section [7] we advocate for a scenario in 
which dust reddened quasars are a short-lived phase in a 
merger-driven co-evolution of galaxies and AGN. In Sec- 
tion[8]we summarize our findings and present ideas for re- 
solving remaining open questions in the future. Any time 
we present the results of cosmological calculations we 
use the following parameters: Hq = 70 km s~^ Mpc""'^, 
Qm = 0.30, and Qa = 0.70. 

2. FINDING RED QUASARS 

In iGlikman et all ()2004 hereafter Paper I) we sought 
to identify the most elusive members of the reddened 
quasar population by searching for objects that were 
detected in the radio and near-infrared, but were un- 
detected in a flux-limited optical survey. We matched 
the 2000 July Faint Images of the R adio Sky at Twenty 
Centimeters (FIRST) radio catalog (|White et al.l 119971) 
with the Second Incremental Point-Sour ce catalog from 
the T wo Micron AU Sky Survey (2MASS: ISkrutskie et all 
|2006( ) and selected objects that lacked an optical coun- 
terpart in the Automatic Plate Measuring (APM^^ ma- 
chine scans of the first generation Palomar Observatory 
Sky Survey (POSS-1) plates. Figure [T] shows the 69 can- 
didates selected in this effort; 54 objects are spectroscopi- 
cally identifiecQ of which seventeen are red quasars. Ap- 
proximately 50% of the objects defined by R — K > 4 
and J — K > 1.7 are heavily obscured quasars. The other 
~ 50% are galaxies exhibiting some form of activity such 
as a starburst or narrow-line AGN activity. Objects with 
bluer J—K colors tend to be late-type stars; sources with 
R — K < 4 tend to be elliptical galaxies due to the large 
4OOOA break, but relatively flat spectral energy distribu- 
tion (SED) redward. 

To further explore this empirically-determined color 
selection, we compute the colors tracks of reddened 
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Paper iT reported 50 spectroscopic observations, subsequently, 



However, see [Lawrence He Elvis! 1)20101 ) for a detailed discus- 
sion of possible departure from this trend. 



four additional spectra - a galaxy and three M stars - h ave become 
availa ble through the Sloan Digital Sky Survey (SPSS: [York et al.l 
l2000h 



FIRST-2MASS Red Quasars 



3 



quasars in J — X vs. R — K space in the redshift range 
z — 0.1 — 2.5. To do this, we combine the Hubble 
Spac e Telescope ( HST) UV composite quasar template 
from lTelfer et al.l (|2002l l w ith the optical-t o -near-infrared 
composite spectrum from IGlikman et al.l (|2006f ) to pro- 
duce a template spectrum that spans 300A to 3.5^m. 
We convolve the template with R, J and Kg transmis- 
sion curves to produce the solid blue line in Figure [TJ 
This curve is nearly identical to the quasar colors pre- 
sented by (jHewett e t al. 200ff) in a study of the optical- 
through-near infrared colors of astrophysical objects in 
the United Kingdom Infrared Te lescope (UKIRT) In- 
frared Deep Sky Survey fUKIDSS: lLawrence et al.l[2007l) 
photometric system. We then redden the template us- 
ing a Small Magellani c Cloud (SMC) re ddening law 
(jFitzpat rick 1999; Gor don fc Clavtonlll998l and see sec- 
tion 15. ip and convolve the resultant spectra to produce 
the yellow, orange and red solid curves corresponding 
to E{B ~ V) — 0.5,1.0 and 1.5, respectively. We also 
explore the effect of host galaxy contribution on the col- 
ors of reddened quasars by adding an unreddened host 
galaxy template that contributes 50% of the overall con- 
tinuum between 0.9 — l.O/inF^ on top of the reddened 
quasar. The dotted lines show the effect of elliptical host 
galaxy, which are the typical hosts of lo w redshift quasars 
(|Duniop et al|[2003l: iFlovd et al.|[200l c.f.,). We also ex- 
plore the effect of a host with younger stellar popula- 
tions by adding an Sb-type galaxy (das hed line). Both 
compo site spectra are taken from the iMannucci et al.l 
(|2001D optical-through -near-infrared galaxy template li- 
brary. We also explored the effect of 70% host galaxy 
contribution, which we do not plot so as not to busy 
the figure. In all cases, the models stay within our color 
cuts, including the models of quasar plus Sc host galaxy, 
which shift the color tracks toward bluer R — K colors, 
but not by enough to remove the sources from our sample 
selection box. 

Having found an efficient color selection criterion for 
reddened quasars, we expanded our search in the same 
area of sky, but includ ed optic al detections. We ex- 
panded the sample from iPaper II by imposing the afore- 
mentioned color cuts on FIRST-2MASS sources but in- 
cluded objects with o ptical detecti ons in the Guide Star 
Catalog II (GSC-II; iLasker et aLll2008[ ): this selection 
addresses the fraction of quasa rs missed due to optical 
reddening (jClikman et al.ll2007l hereafter Paper II). The 
2000 July FIRST catalog has 2716 deg^ areal overlap 
with the 2nd 2MASS incremental release. In this area 
we found 156 quasar candidates. Spectroscopic identifi- 
cation of 77% of these objects resulted in a sample of 52 
red quasars. We used this sample to study their redden- 
ing properties and estimate the overall fraction of quasars 
that are reddened by dust. Because of the small size of 
that sample, as well as the shallow limit of the 2MASS 
survey compounded by the spectroscopic incompleteness 
at the faintest K magnitudes, we could only estimate 
the fraction of red quasars for objects with unreddened 
K < 14.0. We found that at this bright end, red quasars 




T his assumption is based on the findings from IGlikman et al.l 
112009) that a 2L* host galaxy contributes 6% of the total contin- 
uum of an average quasar at these wavelength. Here we increase 
the contribution to 50% to provide a generous limit on the effect 
of the host. 



Fig. 1. — An updated version of Figure 4 from lPaper II showing 
the colors of FIRST radio sources with 2MASS matches and no 
optical counterpart in the POSS-I survey. We find that the color 
cuts R—K > 4 and J—K > 1.7 select quasars with a 50% efficiency, 
avoiding galaxies, which are bluer in R — K and red dwarf stars 
(e.g., M stars), which are bluer in J — K. We also plot color tracks 
for quasars at 2 = 0.1 — 2.5 with different amounts of reddening 
(solid lines) and with the host galaxy contributing 50% of the flux 
at 0.9-1/im. The reddenings are defined in the lower left legend 
(and see text for details). We also plot the effect of a host galaxy 
on the colors of reddened quasars for different amounts of reddening 
with an elliptical host galaxy (dotted line) and an Sc type galaxy 
(dashed line). Although the presence of a star-forming host shifts 
the colors of reddened quasars blueward in R — K color, especially 
for more heavily reddened nuclei, their colors remain within our 
selection box. 

make up 56 ± 16% of broad-line-emitting quasars. 

Taking advantage of the superior depth and im- 
age quality of the Sloa n Digi tal Sky Survey (SDSS; 
lYork et al.ll2000D . lUrrutia et al.l (2009) developed a red 
quasar sample using FI RST, 2MASS and the fifth data 
release of SDSS (DR5; lAdelman-McCarthv et all 120071) 
with slightly modified color cuts from lPaper III : ~— K > 
5 and J — K > 1.3. Be cause of t he overlap of the areas 
and color-selection with Paper II, 28.7% (35/122) of the 
candidates including 42% (24/57) of the red quasars in 
lUrrutia eT al. (2009) were previously identified spectro- 
scopically in Paper H . 

In this paper we more than double our sample of red 
quasars using the full 9033 deg^ overlap between the 
FIRST and 2MAS S survey s . This sa mple uicludes most 
of the objects in IPaper II . IPaper III and lUrrutia et al.l 
(|2009l) and has better spectroscopic sampling at the faint 
end (based on 2MASS K magnitudes) yielding a catalog 
which includes some of the most luminous - yet previ- 
ously unidentified - quasars in the Universe. 

3. SELECTING THE FIRST-2MASS RED QUASAR 
SAMPLE 

Here we describe in detail the construction of our 
red quasar sample, which depends on three wavelength 
regimes: (1) candidate quasars are radio sources, to avoid 
confusion with stars; (2) we require that our quasars 
be bright at 2 ^m, since we are searching for red ob- 
jects; (3) we use optical data to select objects with red 
optical-to-near-infrared colors. In addition, on spectro- 
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scopic followup, we define a red quasar as an object whose 
spectrum shows at least one broad emission line, with 
V > 1000 km s~^, and whose reddening is measured at 
EiB-V)>0.l. 

The FIRST Radio Survey (|Becker et all [l995l ) has 
mapped 9033 deg^ of the sky at 20 cm with the Very 
Large Array (VLA) in the B-configuration. Coverage in- 
cludes ~ 8400 deg^ in the north Galactic cap and ^ 600 
deg^ in the south Galactic cap. With the antennae in this 
configuration, the survey's resolution is roughly 5" and 
has 0'.'5 positional accuracy. The 3-minute snapshot inte- 
gration time yields a typical rms of 0.15 mJy. The 2003 
April 11 catalog has ~ 810, 000 sources b righter than 
the s urvey's detection threshold of 1 mJy (jWhite et al.l 
[1991 . ^ 

Th e 2MASS near-infrared survey (jSkrutskie et al.l 
|2006( ) used two telescopes, one at Mt. Hopkins, Arizona 
and one at Cerro Tololo, Chile, to image the sky in three 
bands, J (1.24 ^m), H (1.66 /xm), and Ks (2.16 /Ltm), 
between 1997 and 2001. This four-year effort produced 
a Point Source Catalog (PSC) with 4.70 x 10* sources. 
The survey is 99% complete to a lOcr magnitude limit 
of J = 15.8, H = 15.1 K = 14.3. Fainter sources are 
included in the catalog (down to K < 16) but with lower 
signal-to-noise ratio and completeness. The 2MASS sur- 
vey has a resolution of 4" and an astrometric accuracy 
better than 0'.'5. 

The Guide Star Catalog II (GSC-Il lLasker et al.|[200l 
is an all-sky optical catalog produced by scanning the 
second-ge neration photogr aphic Palomar Observatory 
(POSS-II IReid et al.l [1991 and UK Schmidt Sky Sur- 
veys at 1" resolution and generating a catalog with 
positions, magnitudes and morphological classifications. 
Wc employed the CSC 2.2.1 catalog which reaches the 
POSS-II photographic plate limits of F < 20.8(0 and 
J < 22.5(T1. The astrometric accuracy of GSC-II is bet- 
ter than 1". 

We matched the 2003 April version of the FIRST cat- 
alog and the 2MASS All-Sky Point Source catalog using 
a 2" radius. In Paper I we showed that a 2" match be- 
tween FIRST and 2MASS has a 0.3% false detection rate 
and is the radius at which background contamination is 
minimized, while preserving completeness. This criterion 
yielded 66,953 matches, including 66,728 unique FIRST 
sources - a catalog we refer to henceforth as F2M. We 
applied the infrared color cut {J — K > 1.7) to this list, 
which reduced the number of candidates to 8443 sources. 
We matched these sources, using the FIRST positions, 
to the GSC-II catalog with a search radius of 2". There 
were 8242 matches within 2" of the FIRST positions in 
F2M. Of these, 8075 are unique FIRST sources, which 
may have two or more CSC and/or 2MASS matches. 
We kept only the closest match to each FIRST source in 
our candidate list. 

In order to apply the color cuts to the 8075 matches 
we need to take into account sources that have been de- 
tected in one optical band but not in another. We re- 
placed the undetected magnitudes for such objects with 
their respective plate limits, B — 22.50 and R — 20.80. 

Photographic F band is equivalent to a ii magnitude, peaking 
at 6750 A. 

Photographic J band is equivalent to a _B magnitude, peaking 
at 5750 A. 



In cases where bright objects were detected on plates 
observed through the Photographic J bandpass but were 
undetected in Photographic F, artificially red colors may 
have been assigned. We examined finding charts of all 44 
such sources and found that 8 were indeed too faint to be 
discerned on the Digitized Sky Survey images; we keep 
them in the candidate list. We obtained magnitudes for 
the remaining 36 s ources using the SDSS-DR8 catalog 
(|Aihara et al.l [20Tl|) which contributed r magnitudes - 
32 were too blue to remain in our candidate list. This 
procedure found 347 FIRST-2MASS sources with CSC 
II matches within 2" obeying our color criteria. 

In Section [3. II we fully quantify the matching statistics 
between F2M and GSC-II, but to finish this section we 
outline some broad trends. Another 358 F2M sources 
had no match in GSC-II within 2". We matched these 
sources to a larger search radius of 6" to find the clos- 
est counterpart. We examined the Digitized Sky Sur- 
vey (DSS) images of these objects and removed 67 ex- 
tended objects (e.g., saturated stars and large galaxies) 
whose FIRST and GSC-II coordinates were separated by 
more than 6", but which were clearly associated with the 
FIRST source. Thirty-six sources had no match within 
6", showing only blank sky on the DSS i mage; el even 
of these objects were previously reported in [Paper 1 . We 
include these sources in our candidate list, adopting the 
plate limits as lower limits to their B and R magnitudes. 

We examined the DSS images of the 247 objects that 
had GSC II matches within 6" but not within 2", to 
search for objects whose real optical counterparts may 
be fainter than the POSS-II plate limit, but which may 
have had chance coincidence matches to another source 
within 6" of their FIRST positions. Many of these ob- 
jects were large galaxies or blended groups of sources that 
overlapped the radio and infrared positions but whose 
centroids implied optical positions separated by 2" and 
6". We removed these sources and retained the remain- 
ing 94 optical objects which appeared isolated. We de- 
scribe below the matching statistics between F2M and 
GSC-II which we use to determine the reliability that a 
F2M - GSC-II association is real. 

3.1. F2M Matches to GSC-II sources 

To determine the completeness and contamination 
statistics of our samp le, we follow the procedure of 
iMcMahon et al.l (|2002D in determining the chance coin- 
cidence rates as a function of separation between F2M 
and GSC II. We matched our 8443 F2M sources to GSC 
II out to 12", keeping only the nearest match (there were 
8354 unique matches). To model the background coin- 
cidence rate, we shifted the 8443 F2M positions by 5' 
to the north and rematched this list to GSC-II, again 
out to 12". We determine the mean source density, 
Popt = 7.09 X 10"** arcsec"^, and the variance in the den- 
sity, o-Qpt = 1.45 X 10"'' arcsec"'', using the formalism 
outhned in §5 of IMcMahon et al.l (|2002D . Figure [2] shows 
the distribution of match separations plotted in 0"2 bins. 
Integrating under the histogram of real matches and the 
chance coincidence curve to 2" determines that 99.8% of 
these matches are physically associated; there are 8075 
real matches and 14 random associations. This result 
implies that of the 343 objects in our candidate list, 0.59 
(or 1) objects are chance coincidences. Integrating be- 
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Fig. 2. — Histogram of separations for F2M-GSC II matches, 
binned by 0'.'2 and normalized by the annular area of each bin 
and number of F2M sources (8443). Only the closest matches are 
included. The inset shows the chance coincidence rate {dashed line) 
at large radii. 

tween 2" and 12" implies that 279 real and 35 coincidence 
matches would be added, suggesting our completeness is 
95.6% (8075/8443). 

Figure [3] shows the separation distribution for point- 
like radio morphologies (left), and extended radio mor- 
phologies (right). We define "pointlike" to be FIRST 
sources with Speak/ S integrated > 0.9. The insct shows 
the chance coincidence rate (dashed line) compared to 
the matched sources. These figures demonstrate that ex- 
tended radio sources are more likely to have real matches 
at larger radii than point sources. We use these results 
to determine the probability of chance coincidence for 
each of the 94 sources whose nearest GSC II match lies 
between 2" - 6" away. We ascribe a probability to each 
source for being a chance coincidence based on its ra- 
dio morphology and separation from the FIRST posi- 
tion. Adding the probabilities of the first ten objects 
(sorted from lowest to highest probability) in this list sum 
to 891%, or ~ 9 chance coincidences. We assume that 
these objects are unmatched sources, and so we adopt the 
POSS-II plate limit for their B and R magnitudes. 

The final candidate list, therefore, contains 395 
sources, 347 of which have detections in the GSC-II cata- 
log, and 48 of which are fainter than the sensitiv ity limit 
of GSC-II. This list includes 31 objects from iPaper II 
that had no detection in the APM catalog and have 
J — K > 1.7 and R — K > 4. A flowchart of our selec- 
tion process is shown in Figure |4l and the list of quasar 
candidates and their attributes is provided in Table [T] 

3.2. Spectroscopic Observations 

Spectroscopic observations of our red quasar candi- 
dates were performed in the optical and/or near-infrared 
over the timespan of ten years with additional spectra 
from archival and published sources in the literature. In 
total we obtained spectroscopy for 316 red quasar can- 
didates. The nature of these spectra is highly heteroge- 
neous, and we use the spectra for two main purposes (1) 
object identification (i.e., galaxy, quasar) and redshift 
determination; (2) for the qua sars, redd ening determina- 
tion (see §5.1 below and §5 of IPaper H ). 

A total of 209 optical spectra were obtained at the W. 
M. Keck Observat ory with the Echellette Spectrograph 
and Imager (ESI; iSheinis et al.l 12002. ) and the Low Res- 



olution Imaging Spectrograph (LRIS; lOke et aI]|1995D . 
the Palomar observatory's 5-meter Hale telescope with 
the Double Spectrograph, and the Lick Observatory's 3m 
Shane telescope with the Kast spectrograph. Additional 
spectra were taken from the SPS S and SDSS-III archives 
(|Gunn et al.l[T998. .2006: Eisenstein et al.llMll ). includ- 
ing the Baryon Oscillation Spec troscopic Survey quasar 
survey (BOSS: IRoss et al.ll20lll) . Column (14) of Table 
[T] lists the origin of the optical spectroscopy. 

The 170 near- infrared spectra (covering the wavelength 
range ~ 0.9 — 2.5 um) were obtain ed at the NASA 
IRTF with SpeX (Ravner et al."200 3D, the Pa lomar 5- 
meter Hale with TripleSpec (Herter et al.ll200"8l) or with 
TIFKAM on the 2.4 m MDM telescope. We reduced 
the SpeX and TripleSpec data using the Spextool soft- 
ware which is designed for cross-d ispersed near-infrared 
spectroscopy ( Gushing et al."2004). The procedure in- 
cludes flat fielding, sky subtraction, extraction and co- 
addition of individual exposures, and wavelength cali- 
bration. The spectra are also corrected for telluric ab- 
sorption using the spectra of nearby AOV stars that were 
obtained di rectly before or after each target observation 
(jVacca et a l. 2003). Column (15) lists the origin of the 
near-infrared spectroscopy. 

A total of 80 objects have both an optical and a near- 
infrared spectrum. Having both spectral regions is es- 
pecially useful in cases where the nature of an object 
is not readily identifiable. For quasars, a near-infrared 
spectrum enables identification of heavily dust-obscured 
objects whose rest-frame ultraviolet emission lines are 
beyond detection due to extinction. On the other hand, 
the optical spectrum better-constrains our estimates of 
E{B — V), since the rest- frame UV is more sensitive to 
dust than longer wavelength light. For galaxies, a near- 
infrared spectrum has fewer identifiable features than in 
the optical and is often difficult to classify. Thirty-four 
candidates with only a near-infrared spectrum could not 
be classified because they lacked strong emission lines or 
other features. 

3.3. Spectroscopic Identifications 

Of the 316 objects with spectra (out of 395 candi- 
dates) we have been able to identify 263, thus our survey 
is 80% spectroscopically complete (67% when consider- 
ing only identified spectra). The distribution of spec- 
troscopic completeness as a function of K magnitude is 
shown in Figure |5l The distribution of all F2M candi- 
dates is shown in the unshaded histogram. Overplot- 
ted are histograms of spectroscopically observed objects 
(forward hash), objects with spectra whose type is iden- 
tified (cross hash), and quasars (black) found in our sur- 
vey. We are 96.6% spectroscopically complete for candi- 
dates brighter than K = 14.75; our completeness drops 
to 72.8% for K > 14.75. 

There are 135 spectroscopically confirmed quasars in 
our survey, defined as having at least one broad emission 
line, with a width of u > 1000 km s~^ based on a Gaus- 
sian fit. We use the emission lines to assign a redshift 
to each quasar. Three additional candidates are identi- 
fied as quasars in the literature; we list their references 
in Table |TJ We present the spectra of 120 red quasars, 
defined as having E{B - y) > O.lQ in Figure |6| in order 

The determination of E(B — V) is outlined in Section IS.ll of 
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Fig. 4. — Schematic diagram of our selection process. 



of decreasing redshift. 

As noted above, as well as in iPaper J . iPaper III and, 
e.g., iRawlings et al.l (|1995l ). reddened objects may only 
reveal narrow lines in their optical spectra (which show 
rest-frame UV emission, depending on the redshift), 
while a near-infrared spectrum may be needed to re- 
veal broad emission lines confirming the presence of a 
quasar. For example, nearly all of the red quasars in 
lUrrutia et al.l (2009) were identified from optical spec- 
tra. F2MJ104043.66+593409.55 is listed as a narrow-line 
AGN (NLAGN) at z = 0.147 in Table 1 of lUrrutia et"all 

this paper as well as §5 of lPaper . 




Fig. 5. — Histogram of the K magnitudes for the F2M candidates 
{unshaded) . All spectroscopically observed objects are overplotted 
in the forward-hashed areas. Objects whose spectrum yielded an 
identification are overplotted in the cross-hashed histogram, and 
the quasars are overplotted in black. Our overall spectroscopic 
completeness is 80%, but rises to 96.6% for K < 14.75. 



(|2009[) : however, our near-infrared spectroscopy of this 
object reveals strong, broad Paschen a 18756A, Pa/3 
12822A, and Pa7 1094lA atop a very red continuum, 
confirming that it is a Type 1 quasar. In fact, of the 
eleven objects lab eled NLAGN b ased only on optical 
spectroscopy in Ur rutia et al.l (12009). six are re-classified 
as quasars in this work based on near-infrared spec- 
troscopy. 

This strongly suggests that many of the objects classi- 
fied here as NLAGN from only optical spectra will reveal 
broad emission in the near-infrared and argues for near- 
infrared spectroscopy of all NLAGN with J — K > 1.7 to 
recover any missed red quasars and to better determine 
the fraction of quasars that show no broad lines in the 
optica l. On the oth e r han d, there are five red quasars 
in the lUrrutia et al.l ()2009D sample that lie outside our 
color cuts, with 1.3 < J — K < 1.7. This amounts to 
~ 9% of that sample (5/56), which, if extended to our 
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Fig. 6. — F2M candidates identified as red quasars, having at 
least one broad emission fine and E{B — V) > 0.1, ordered by red- 
shift. The dotted fines show expected positions (in ^tm) of promi- 
nent emission fines in the optical and near-infrared: Lyo 1216, 
N V 1240, Si IV 1400, C IV 1550, C III] 1909, Mg II 2800, 
[O II] 3727, H(5 4102, H7 4341, H/3 4862, [O III] 4959, [O III] 5007, 
Ha 6563, Pa7 10941, Pa/3 12822 and Paa 18756A. The fuU atlas is 
presented at the end of the paper. 



sample, suggests that we are missing ~ 11 red quasars 
because of our J — K color cut. 

We also find nine BL Lac objects: highly variable, 
featureless blue objects with bright radio flux densities 
(>_100 mJy). Two of these AGNs are repeats from 
iPaper IB . The remaining seven have been identified as 
BL Lacs elsewhere; references are li sted in co lumn (16) 
of Table [H As demonstrated in §8 of lPaper^ . these blue 
objects are selected in a survey targeting red objects be- 
cause of significant variability between the epochs of opti- 
cal and infrared observations, where the former detected 
the source in a faint state and the latter was viewed with 
the object in a bright state. We exclude these objects in 
our analysis for the remainder of this paper. 

Fifty-two objects reveal narrow line emission in their 
spectra and are identified as either AGN and/or star- 
forming galaxies. For objects with an optical spectrum, 
wc model the H^, [O III], [O I], [N II], Ha, and [S II] 
lines with Gaussian profiles and pl ot their line ratios on 
Baldwin-Phillips-Terlevich (BPT; [Baldwin et all [TOSl 
line-ratio diagrams to determine the source of the ion- 
izing flux (AGN, star formation, composite object show- 
ing contribution from both AGN and star formation) 
for each spe ctrum (Fi gure [71). This includes seventeen 
objects from [Paper IB . Three narrow-line emitters, e.g., 
F2MJ104902.95-h401031.6 at z = 0.715, only have spec- 
troscopic coverage of H/3 and [O III] from an optical spec- 
trum. With log([0 III]/H/3)=0.75, the nature of this ob- 
ject's emission is not clear, but it is likely to have some 
contribution from star formation and is probably a com- 
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Fig. 7. — Line diagnostics for the n arrow-line-emitting sp ectra, 
plotted based on the classifications of IKewlev et al.1 II2006I) . The 
pure-star formation boundary is shown with a dashed fine in the 
leftmost panel, while the Seyfert-LINER divide is shown with a 
dashed line in the rightmost panel; LINERs have lower [O IIl]/H/3 
values than Seyferts. The left panel shows the [O IIl]/H/3 vs. 
[N II] /Ha diagnostic. In this panel, objects below and to the left 
of the dotted line are pure starbursts. The middle panel shows 
the [O IIl]/H/3 vs. [S II] /Ha diagnostic. The right panel shows 
the [O III]/H^ vs. [O I]/Ha diagnostic. Objects to the left and 
below the solid line are starburst dominated (filled circles), while 
objects above and to the right of the line are AGN dominated (open 
circles) . 

posite object. We list the line diagnostics from our mea- 
surements in Table [H including, for completeness, the 
sources already analyzed in Paper II. 

Thirteen of the narrow-line-emitting objects have only 
a near-infrared spectrum. We cannot distinguish be- 
tween their source of ionization based only on narrow 
Paschen lines, for which there are no BPT-diagram-like 
diagnostics. However, the presence of bright radio emis- 
sion suggests that they may host AGN. 

Fifty-six of our objects are galaxies with no obvious 
emission lines. We determine their redshifts by cross cor- 
relating the ir spectra with ellipt ical and SO galaxy tem- 
pl ates from [Kinney et al.[ ()19960 . following the method 
of [Tonrv fc DavisI (|1979l) 

There are five stars in our sample, which is a fraction 
consistent with the previous F21VI red quasar samples. 
We assign stars a redshift of in Table [1] 

4. THE OBSERVED PROPERTIES OF RED 
QUASARS 

4.1. Comparison Samples 

Our final F2JVI red quasar sample consists of 120 ob- 
jects spanning z ~ 0.1 — 3. In contrast to optically- 
selected quasars, these objects have optical-to-near- 
infrared SEDs that rise toward longer wavelengths. Are 
we simply finding the tail end of an intrinsic distribu- 
tion of quasar slopes (e.g., the red quasars discussed in 
Rjxiiards et al. 2003) or axisynimetrically dis tributed nu- 
clear dust (e.g., type 1.5 quasars presented in lSmith et al.l 
[2002f ) or do these heavily reddened quasars represent a 
separate population of objects, possibly associated with 
a phase in quasar e volution (i.e., the "blowout" phase 
[Hopkins et al. I [20081 )? To address this question, we must 
compare the F2M quasars with an appropriate optically- 
selected radio-detected quasar sample. 

The most natural comparison set for the F21VI red 
quasars is the optically-selected FIRST Bright Quasar 
burvey CFBQS: [Gregg et al.lfT996l) . This survey searched 
for quasars in matches between FIRST radio sources and 
sources in the APM scans of the first generation Palomar 
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all sky survey (POSS-I). The main survey contains 636 
quasars over 2682 deg^ and imposed an £' = 17.8 magni- 
tude limit on the APM sources as well as an O — i? < 2 
color criterion (FBQS II; White et al. 2000). A deeper 
segment of the survey covered 589 deg^ and found 321 
quasars wit h the same color c ut, but with E < 18.9 
(FBQS III; IBecker et all [200l . Isolating the FBQS 
quasars that are detected in 2MASS forms an optically- 
selected sample of quasars that exist in the A'-band flux- 
limited survey; there are 503 and 134 2MASS matches to 
FBQS II and III, respectively. 

We also construct a, subs ample of quasars from the 
SDSS (|Schneider et al.l[2010l) that are detected in both 
FIRST and 2MASS to compare the efficiency of the SDSS 
quasar selection, which is more sensitive to reddened ob- 
jects than simple UVX color selection because it includes 
FIRST-ta rgeted sources as well as outliers from the stel- 
lar locus (jRichards et al.|[2"002[ ). We estimate the effec- 
tive area for radio-detected SDSS quasars to be the area 
of the 2003 April 11 FIRST catalog release (9033 deg^) 
minus ~ 300 deg^ which appear to be lacking SDSS cov- 
erag43- We then apply the 7 7.4% spectroscopic com- 
pleteness (jRichards et al.|[2006[) to this area for an effec- 
tive area of 6770 deg^. 
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Fig. 8. — Redshift distribution of as a function of B — K color 
of F2M quasars {red triangles) compared with FBQS II (blue open 
circles) and III {filled circles). The F2M and FBQS surveys reach 
the same redshift ranges, making them good comparison samples. 



4.2. Observed Surface Density 

Figure [8] plots B — K color versus redshift for F2M 
and FBQS quasars and shows that the quasars found in 
the F2M survey span the same redshift range as FBQS 
quasars. In Figure [HI we plot a histogram showing the 
distribution of the B - K color for FBQS II and III 
(solid and dashed line histograms, respectively) next to 
the same for the F2M quasars (shaded histogram). The 
histograms have been normalized t o each survey ' s cov - 
era ge area fo r ease of comparison. Webs ter et al.l ()1995[ ) 
and lPaper II showed that optically-selected quasars have 
typical B — K ~ 2.5 color, which is not a strong function 
of redshift, as can be seen in Figure |S1 Therefore B ~ K 
color can be used as a rough proxy for reddening. 

Figure [9] shows the distribution of the B — K colors of 
FBQS and F2M quasars. The overlap between the two 
samples is small, by construction, but it appears that 
the red quasar color distribution (shown in more detail 
in the inset) does not simply fill in the tail end of the op- 
tically selected B — K color distribution. This would be 
expected if there were no reddening in quasars; the colors 
would then represent the intrinsic variation of their op- 
tic al continuum sl opes. A similar observation was made 
bv [Richards et all |2003) who examined the distribution 
of the relative coloiFl o f SDSS quasars. In their Figure 
3, [Richards et al.l (|2003D note that there is an excess of 
sources with red A(g — i) colors compared to a Gaussian 
fit to the distribution defined by its peak and blue wing. 
They conclude that that there are some quasars that are 
intrinsically red because their continuum has a flatter 
spectral index, but exist within the normal distribution 
of spectral indices for quasars. The quasars far in the 

We estimate this area by computing the spherical area of this 
lune from its four corners. 

Defined as the residual color of a quasar at a given redshift 
after sub tracting the median quasar color at the same redshift 
IIHicharcIsil2001]) . e.g., A((; — i) and correcting for Galactic extinc- 
tion. 
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Fig. 9. — B — K distribution of F2M quasars compared with 
FBQS II and III; arrows denote bins which are lower limits. The 
distribution of the optically selected (FBQS) samples peak near 
B - X ~ 2.5, though the distribution of FBQS-III (dsahed line), 
which goes 1.1 magnitudes deeper than FBQS-II, is slightly shifted 
redward demonstrating that a fainter optical flux limit will picks up 
some redder sources. The F2M quasars appear to have a separate 
distribution, and not just an extension of the red tail of FBQS 
quasars, suggesting that these objects are not merely a tail end of 
the intrinsic distribution of quasar colors. 



tail of the /^{g* —i*), however, are reddened by dust. We 
address this in the context of red quasars representing a 
phase of quasar evolution in SjT] 

Figure [TU] shows the observed number counts of the 
F2M quasars (red squares and line) compared with the 
number counts of FBQS quasars (plus sign for FBQS 
II and "x" symbol for FBQS III connected by a dot- 
ted line) . Table |3] lists the number counts for the F2M 
quasars plotted in Figure 1101 To determine the fraction 
of quasars that are red, we integrate the F2M curve and 
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Fig. 10. — Spatial density of quasars on the sky of F2M red 
quasars (filled red squared) compared with FBQS II and III (dot- 
ted lines) and S DSS quasars with FIRST detections (dashed Hne; 
[Schneider et al.]|20 10). Radio plus optical selection for quasars 
misses 11 it 2% of quasars because of reddening. The open violet 
squares show the spatial density F2M red quasars after correcting 
for X-band absorption. Beyond K = 14.5 incompleteness from 
the 2MASS flux hmit begins to become apparent. Comparing the 
space distributions of F2M quasars with FBQS we find that red 
quasars make up 21 ± 2% of radio-selected quasars with K < 14.5. 



an average of the FBQS curves out to K 
compute their ratios: 



N{K < 15.5)f2M 
N{K < 15.5)fbqs 



X 100. 



15.5 and 



(1) 



Comparing the total areal densities, we find that FBQS 
quasars miss 11 ± 2% of radio-selected quasars because 
of their red colors. 

We also plot the number counts of the radio-detected 
SDSS quasars described above with a dashed line. Since 
many F2M quasars had SDSS spectroscopy, there is con- 
siderable overlap between the samples. Therefore, the 
fraction of red quasars missed by SDSS is smaller. De- 
spite this, Urrutia et al. (2009) showed that ^ 60% of 
the red quasars in their sample are not selected for 
the m ain z < 3 sample (e.g., see IRichards et all 120021 
I2004D . Furthermore, the completeness of SDSS drops 
significantly in redshift ranges that correspond to colors 
overlapping the stellar locus, which is most pronounced 
for 2.7 < z < 3.2 - at the edge of our survey - and 
1-5 z < 3.2 for Broad Absorption Line (BAL) quasars 
with large absorption troughs (we elaborate more on this 
in Section in]). For the F2M quasars in this work, we find 
that 89 of our candidates have SDSS spectraPI 

It is the targeting of FIRST matches in SDSS that 
increas es the probabil i ty of a spectroscopic selection in 
SDSS (jRichards et al.l 120021) . This feature means that 
there is a radio-dependent bias of red quasar spectra 
in the SDSS spectroscopic database. Radio-quiet red 
quasars (which fall below the FIRST detection thresh- 

Note that we hst SDSS as the main source of spectroscopy 
for only 61 candidates in Table [T] because we only list SDSS as 
the source of the optical spectroscopy in the absence of superior 
quality spectra from our own observations. 



old) would therefore be less likely to be spectroscopically 
targeted by SDSS, resulting in a large missed fraction 
of red quasars in the analyses of quasar properties (e.g., 
luminosity functions, clustering) determined from this, 
and other optically-selected samples. 

4.3. Radio Properties 

Since the quasars in this survey are selected in the ra- 
dio, we examine their radio properties to determine if 
radio-selection introduces a bias. Although quasars were 
originall y discovered as optical counterparts to radio- 
sources (jSchmidtl 119631 ) it was not long before it was 
realized that only ^ 10% of quasars are strong ra- 
dio emitters (jSandagel I1965D . The radio properties of 
quasars are often described using the radio-loudness pa- 
rameter defined as the r atio of radio to op tical power 
{R = ffl.4GHzyf(B^.t): [Stockeet al.l[T99a ). The tra- 
ditional definition of a "radio-loud" quasar is one hav- 
ing ^ R* > 10, where R* is the if-corrected radio- 
loudness parameter. The majority of FBQS quasars exist 
in 3 < R S. 100 an d span a "radio intermediate" regime 
(jWhite et al.1 120001) . 

In PapcrJ we computed the radio-loudness parame- 
ter for the seventeen F2M quasars in that sample using 
the extinction corrected i?-band magnitude, Bint- We 
found that these sources also lay in the radio intermedi- 
ate regime, with only two quasars having R > 100. 

Here we examine the radio-loudness for the full F2M 
quasar sample. We use a more current defini tion of 
the radio-loudness parameter, as defined bv Jvezic et al.l 
11002), which uses the difference between the mag- 
nitude equivalent of the FIRST flux density (t = 
-2.51og(Fint/3631Jy)) and the SDSS magnitude in a 
given bandpass, m, both on the AB magnitude system: 



Rm = l0g(f;adio/^'optical) = 0.4(r71 - t) . 



(2) 



We use the g-band, which is the closest to the previously- 
used J5-band, for this calculation. This definition of Rm 
does not include a fsT-correction, which avoids adding 
uncertainties introduced by the unknown radio spec- 
tral indices of these FIRST-selected sources, as well as 
the already-uncertain reddening corrections (see Section 
(51) which would be compounded by the uncertainty of 
adding a JC-correct ion t o the optical mag nitude (see 
iGlikman et al.l 120111 and iCroom et al.l 120091 for discus- 
sions of the uncertainties introduced by iC-corrections in 
quasars) . 

Figure [11] left, shows the distribution of Rg for the 
F2M red quasars compared with FBQS quasars. The 
dashed line shows Rg for the F2M red quasars using their 
apparent g-band magnitudes (uncorrected for reddening) 
from SDSS while the dot-dashed and dot-dot-dot-dashed 
hues arc from FBQS-II and FBQS-III. The uncorrected 
F2M distribution suggests that red quasars are skewed 
toward radio-loud sources. However this is because dust 
extinction suppresses the denominator in Equation [2] ar- 
tificially enhancing Rg. Once we correct for extinction 
(our methodology for dereddening is described in Section 
[S]) the Rg distribution for the F2M quasars (solid line) 
has the same distribution as the FBQS quasars (dash-dot 
line). We see that for both the FBQS and F2M samples 
many are bona-fide "radio-quiet" quasars with Rg < 1. 
Still, most quasars lie below the FIRST detection thresh- 
old and they are even "quieter" in the radio; only 16% of 
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the q uasars in the SDSS quasar catalog (jSchneider et alj 
[MOI overlap the FIRST survey area and are detected 
in F IRST (an d have 2MASS detections). 

In iPaper III we reported contemporaneous 3.6 cm and 
20 cm flux density measurements with the VLA for 44 
F2M quasars. The distribution of spectral indices for 
F2M quasars showed that there are few flat spectrum 
sources (q!i.4ghz/8.3GHz > —0.5, where 5^ oc u"') and 
that the fraction of steep-spectrum sources rises with de- 
creasing 20 cm flux density. The dearth of flat radio 
spectral indices in the F2M quasar sample suggests that 
their red colors are not caused by synchrotron contam- 
ination. This is in contrast to the red quasars found 
in the Par kes Half-Jansky Flat-Sp ectrum quasar sample 
(PHFS iDrinkwater et al.lll997D where iWhiting et al.l 
(POOl) do find a red synchrotron component contribut- 
ing to their red colors. To see if very bright radio sources 
contain some red synchrotron emission that may mas- 
querade as dust-reddening, we plot E{B — V) vs. 20 cm 
FIRST flux density for all F2M quasars in Figure [IH 
right. We see that below ^ 20 mJy there is no cor- 
relation with reddening and that the most radio-bright 
sources tend to be the least reddened. In particular, 
seven of the eight F2M quasars with i^20 cm > 500 mJy 
have E{B — V)< 0.4 suggesting that a red synchrotron 
enhancement is not important in these sources. We ar- 
gue, therefore, that the red colors of F2M quasars are the 
result of dust extinction. The range of spectral indices in 
the F2M quasars also suggests that we are viewing the 
F2M quasars at a variety of orientation angles. 

If the dust that reddens the quasars is located in 
their host galaxies as part of a merger-driven picture for 
quasar/galaxy co-evolution (see Section[7]) then we would 
not expect the radio emission from the quasar to be re- 
lated to dust located parsecs to kiloparsecs away from 
the black hole. Therefore, our findings here should ex- 
tend to the overall quasar population, regardless of their 
radio properties. 

Recent results, however, suggest a possible correla- 
tion between redde ning in quasars and radio emission. 
iWhite et al.l ([2007') found, in radio stacking of quasars 
(the overwhelming majority of which are radio-quiet) 
that redder quasars have high e r med ian radio fluxes. 
In addition, iGeorgakakis et all (|2009[ ) studied the ten 
brightest red quasars in 2MASS, selected without any 
radio constraint, and found that 60% were detected in 
the radio (compared with ~ 10% for all quasars). They 
also reported evidence for high levels of star forma- 
tion based on Spitzer photometry su ggesting that these 
QSOs are young. Follow-up work bv IGeorgakakis et al.l 
(|2012D find that such bright red quasars have spec- 
tral properties in the radio indicative of young radio 
jets. A radi o- independent selec tion of red quasars (e.g., 
with WISE; iWright et al.ll2010l ) is needed to determine 
whether we can extend the conclusions drawn from this 
radio-selected sample to the entire quasar population. 

4.4. Morphological Properties 

In addition to their blue, UV excess SEDs, unobscurcd 
quasars outshine their host galaxies and have a stellar 

22 The PHFS radio flux limit is 500 times higher than FIRST, so 
the contribution of red synchrotron emission in the optical is only 
an issue at high radio luminosities. 



morphology. Many quasar samples have exploited this 
feature, requiring that the images of candid a,te quasars be 
unres olved point sources (e.g., FBQS, 2QZ, ICroom et al.l 
2001: 'Glikman et al.ll2010h . Reddened quasars, however, 
can appear extended in their rest-frame UV and optical 
images as the quasar emission is attenuated relative to 
their host galaxies. 

The SDSS assigns a morphology for all sources de- 
tected in its five imaging filters (type_u, type_g, 
typej:, type_i ,type^). A global morphology, (type), 
which is based on the combined flux from the five pho- 
tometric bands, is also assig ned (iStoughton et al.ll2'00l . 
Figure [T2| shows the distribution of F2M quasar mor- 
phologies, shown as a fraction of the total number of 
sources (shaded bar). For comparison, we also plot the 
SDSS morphologies of the FBQS quasars, which were se- 
lected to have stellar morphologies in the digitized POSS- 
I data (black bar). We see that roughly half of the F2M 
quasars have an extended morphology, implying that sur- 
veys for red quasars that apply a morphological restric- 
tion, especially in the optical, can miss ^ 50% of the red 
quasars. 

We also checked for a color dependence on morphol- 
ogy by examining the classification of F2M quasars in 
the individual filters. We find that the ratio of 'STAR' 
and 'GALAXY' morphologies are distributed similar to 
the global distribution in the u, g, r and i filters. In 
the z band, however, ^ 60% of F2M quasars have a stel- 
lar morphology, consistent with dust-reddening in these 
quasars. 

5. REDDENING IN QUASARS 
5.1. Estimates of Reddening in Quasars 

In IPaper ij we analyzed the reddening properties of 
56 quasars, a subset of the sample in this paper obey- 
ing identical selection criteria. We applied a stan- 
dard extinction law to an optical- t o-nea r- infrared quasar 
template (/o; iBrotherton et al] 12001': 'Glikman et alj 
2006 ) using an SMC red dening law (Fitzpatrick 199^ 
Gord on &: ClavtonI '1998), which has been shown to 
fit reddening in quasars more effectively than a Large 
Magellenic Clo u d (LMC) or Milky Way red dening law 
(jHopkins et al.l l200l iRichards et all f2003l ). i gnoring 
emission lines and noisy regions. Here we follow this 
method for the full F2M quasar sample and use the re- 
sultant E{B — V) values throughout the paper. Since 
we obtaine d additio nal spectra for some of the quasars 
reported in IPaper III , adding in a missing optical or near- 
infrared spectrum, we re-compute E{B — V) for those 
quasars along with the expanded sample. 

In Paper II, we noted several caveats to consider when 
determining reddening: (1) The uncertainty in the scal- 
ing between the optical and near-infrared spectra can 
result in poorly constrained fits, and likely erroneous 
E{B - V) values (e.g.. Figure 10 of IPaper 11 . The scal- 
ing is performed in the overlapping ends of each spec- 
trum's wavelength range, which are the noisiest and 
worst-calibrated. (2) Reddening determined from fits to 
the combined spectrum sometimes differed significantly 
from fits to the optical or near-infrared spectrum of the 
same object. In such cases we typically rely on the fit 
to the optical spectrum, since it is always sampling the 
shortest wavelength light, which is most sensitive to ex- 
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Fig. 11. — Left - Radio loudness distribution for F2M quasars based on observed SDSS g-band magnitudes (dashed line) compared with 
FBQS quasars (dot-dahsed line), as defined in Equation[2l Both distributions have been normalized by their respective survey areas. Based 
on the observed optical magnitudes, which are dimmed because of extinction, the F2M quasars appear to be radio "louder" than the FBQS. 
Once we correct the g-band for extinction, however, the radio- loudness for F2M quasars decreases and appears to be distributed similarly 
to the FBQS quasars. Right- A plot of reddening E{B — y)vs. FIRST 20 cm flux density shows that the most radio bright quasars (filled 
circles have F20 cm > 500 mjy) are the least reddened. This suggests that a red synchrotron component is unlikely to be the cause of the 
red colors of F2M quasars. 
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Fig. 12. — Morphological distribution of F2M quasars based on 
the global SDSS morphological classification (shaded bar). Half of 
these reddened quasars appear extended, making them elusive in 
surveys that select for point sources. For comparison, we show the 
distribution of the same for the FBQS quasars (black bar). Both 
distributions are normalized by the total number of quasars in that 
sample. 

tinction, and therefore best constrains the absorption pa- 
rameters. (3) The presence of host galaxy hght, espe- 
cially in low-redshift, low-luminosity sources, where the 
4OOOA break is occasionally visible, can affect our red- 
dening estimates. We checked for this issue in IPaper IJ 
and found this effect to be negligible. Furthermore, re- 
moving the galaxy in some sources but not others adds 
an inconsistency to our analysis. After all, the galaxy is 
there in all cases, but our ability to quantify its contribu- 



tion to the light depends on many factors, e.g., signal-to- 
noise ratio of our spectrum around the rest-frame 4000A 
break, the wavelength range of a given spectrum, etc. 
We therefore do not remove the galaxy in our reddening 
analysis. 

In HST images of a subsample of F2M quasars 
lUrrutia et al.l (j2008D performed detailed morphological 
analysis and point-spread function (PSF) subtraction for 
these sources and found that the g — I colors of the point 
sources were even redder than their total aperture col- 
ors. The host galaxies also showed evidence for merging 
and interactions, which suggests that star-formation in 
the host galaxies of these sources added blue light to the 
lower-resolution observations (this is corroborated by the 
color-tracks shown in Figure [1] we discuss the implica- 
tions of these HST observations in Section [7]) . We may 
therefore be underestimating of E{B — V) in our quasars. 
Table [3] lists the extinction parameters for 131 quasars 
in our sample that have sp ectra. 

We also investigated, in IPaper I j . the correlation be- 
tween reddening derived from Balmer decrements (i.e., 
the ratio of Ha to H/3) measured from our quasar spectra. 
We only conduct this analysis in cases where both lines 
appear in a single spectrum, to avoid uncertainties intro- 
duced from scaling. We found that E{B — V) measured 
from Balmer decrements was consistent with E{B — V) 
measured from continuum fitting, on average, but var- 
ied significantly from object to object (with a scatter of 
~ 0.5 mag). 

Although Balmer decrements have been used histor- 
ically to measure redde ning along the line o f sight to 
AGN and quasars (e.g., iMaiolino et"alll2001[ ) we argue 
that using the full continuum to model the extinction to 
a source is a far more reliable approach. First, in both 
approaches, E{B — V) is effectively derived by determin- 
ing the relative extinction from a fiux ratio at different 
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wavelengths compared to an intrinsic flux ratio, 



E{B - V) 



1.086 



log 



/(A) 
/o(A) 



(3) 



The same dust law, fc(A), is used for measurements of 
Balmer decrements, replacing fc(A) with k(Ha) — k(H/?), 
/(A) with the measured ratio of Ha/H/3 line fluxes and 
/o(A) with an intrinsic ratio of Ha/H/3 in Equation [3l 
Therefore, in both methods the same formulation is used 
to determine E{B — V) except that more wavelengths are 
sampled when the full continuum is used, which makes 
it a more robust measurement. 

Secondly, both methods require an assumption about 
the intrinsic flux ratios of parts of a quasar's spectrum as 
a function of wavelength (/o(A) versus Ha). The intrin- 
sic shape of quasars' continua are known to have intrinsic 
variation that f ollows a Gaussian distribution (/^ oc 
with aa = 0.30: lVanden Berk et al.]|2001l:lRichards et alj 
[200l . This means that an intrinsically redder-than- 
average quasar may be over-corrected in a continuum 
fit. We partially mitigate this issue by using the FBQS 
quasar composite spectrum (Brotherton et al. 2001) as 
the rest-frame UV to optical portion of the template that 
we use to fit our quasars. Since this composite is derived 
from the same quasar sample to which we are comparing 
our F2M quasars, our analysis remains i n terna lly consis- 
tent. In addition, the IBrotherton et al] (|2001f ) compos- 
ite was shown to be slightly redder (a = —0.46) than 
quasar tem plates derived from optically-selected quasars 
(a — — 0.32: lFrancis et al.lll99l[ ) making it less likely that 
we are overes timating the reddening in our quasars. Fur- 
thermore, in iPaper Ij we explored the effect of varying 
the intrinsic slope of the quasar composite template to 
the ±1(Tq, on the values of E{B — V) derived from contin- 
uum fits; the variation is < 0.1 magnitudes in E{B — V) 
(see §5.3 and Figure 14 of that paper). 

On the other hand, the variation in the intrinsic 
Balmer decrement of the broad line region can have a 
wide range of values spanning at least a factor of two. 
This large range is a result of collisional excitation of 
the Balmer lines and radiative transfer effects at the 
high densities (A^e ^ 10^~^° cm~'^) of the broad line re- 
gion (|Netzeiill975l: iRees et al.lll989[) . The traditionally- 
used intrinsic Ha/ H/3 ratio, derived from Case B recom- 
bination, is 2.88 ('Oste rbroclj 119890 . However, Case B 
recombination is likely ruled out in the broad line re- 
gion, as show n by the Paschen line ratios in the near - 
infrared (e.s.. iGlikman et al.|[200l lOvabu et "all [20091 ) . 
We measure a Balmer decrement for the FBQS tem- 
plate spectrum of 4.091 for the total line profile. With 
such high variations in the intrinsic Balmer decrements 
of quasars, it is not surprising that there is a large scat- 
ter when comparing Balmer decrement-derived redden- 
ings to continuum-based reddening measurements. We 
therefore consider the Balmer decrement to be a far less 
reliable reddening indicator compared with fitting a red- 
dened template to a full spectrum. 

We also explored whether the SMC dust law is the ap- 
propriate reddening law to use for the F2M quasars. To 
do this we fit our spectra with an a verage Large Magel- 
lenic Cloud (LMC) dust law from (iMisselt et al.lll999D 
as well as the Milky Way dust law of ICardelli et al.l 
(fl98a) (CCM, henceforth) updated in the near-UV by 
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Fig. 13. — Comparison oi E(B — V) values derived from different 
dust laws applied to reddened template fits to the F2M quasars. 
We plot the difference between th e SMC and the LMC, Milky Way 
(CCM) and lCalzetti et al.l iwM ) starburst dust laws [AE{B - V)] 
as a function of redshift. We see that the discrepancy increases 
toward higher redshifts where more rest-frame UV light - where 
the dust laws begin to diverge - is used in the fitting. At z > 1 
the discrepancy between the derived extinctions always measure 
the lowest E{B - V) for the SMC dust law. 



lO'Donnelll (|1994l ). In addition, we fit the iCalzetti et al.l 
(1199 4!) dust law for starburst galaxies to explore its 
effectiveness at fitting the continuum of red quasars. 
IGordon et all |2003) performed a comparison between 
the SMC, LMC and Milky Way dust laws and found that 
at rest frame near- UV (- 2500A) through near infrared 
('^ 1/im) the three dust laws are extremely similar. The 
extinction curves begin to diverge shortward of ^ 2000A. 
Figure [T51 shows the difference between E{B — V) derived 
using an SMC dust law and the other dust laws, as a func- 
tion of redshift. We find excellent agreement between 
the values oi E{B - V) derived using the SMC, LMC 
and CCM for most objects, with the standard deviation 
between E{B — V) derived using the SMC dust law and 
the LMC and CCM laws are aE(B-v) = 0.055 and 0.089 
ma gnitudes, respectively . The extinctions measured with 
the lCalzetti et all ()1994f ) dust law have much larger scat- 
ter and poorer agreement with a crE{B-v) = 0.189. We 
see that the discrepancy between the dust laws increases 
at increasing redshifts since we are sampling more of the 
rest-frame UV part of the quasars' spectra. Th is is con- 
sistent with the findings of IGordon et al.l (|2003f ). 

We note that in cases where there is a large discrepancy 
between these two dust laws, the SMC dust law yields 
a lower value of E{B — V). This means that the SMC 
law is a conservative choice when estimating reddening in 
these quasars. We plot in Figure [HI example fits to F2M 
quasar spectra showing instances where the four dust- 
laws produce good fits and excellent agreement (top three 
spectra) and cases where the SMC dust law produces 
superior fits to the others, especially around rest-frame 
A = 2 175 A (bottom three spectra). 

In addition, quasars can sometimes be reddened 
by int ervening absorbers with CO, H2, and HD 
([Noterdaeme et al.,.2010, ) as well as by dust in interven- 
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Fig. 14. — Comparison fits of a quasar template re ddened by four 
dust laws, SMC, LMC, Milky Way (CCM) and the lCalzetti et all 
([1994) law for starburst galaxies to F2M quasars. The top three 
panels show cases where all four dust laws produce similar satis- 
factory fits, and E{B — V) values that agree to within ~ 0.1 mag- 
nitudes. The bottom three panels show that in cases where the 
discrepancy between the fits is large, the SMC dust law produces 
the most consistent fits to the spectra, especially at high redshifts 
where the 2175A "bump" characteristic of the SMC and LMC dust 
laws is not seen. 

ing Mg II absorption systems, whose averag e extinction 
curve s are very similar to an SMC dust law ()York et alj 
|2006|) . iSrianand et~all ()2008D find an LMC extinction law 
in Mg II absorbers in two quasars with E{B ~V)> 0.3. 
Among the F2M sources, we find quasars with narrow 



absorption-line systems whose Zabs 



< 



Absorption 



Hnes that are within 5000 km s of the emission lines 
are considered eit her 'intrinsic' or 'associated' absorbers 
(iFoltz et alJll"986l) . In some of the objects the redshifts of 
the emission lines and the absorption lines are identical 
within their respective errors. Recently, ,Shen fc Menard 
(|2011h have demonstrated that quasars with associated 
absorbers have signatures of enhanced star formation. 
They suggest that these absorbers are possible large-scale 
outflows indicative of a phase in a merger-driven evolu- 
tionary scenario for quasars (see Section [7]). We will ex- 
plore the sub-population of F2M red quasars that show 
broad, associated absorption in their spectra in Section 

El 

5.2. Trends With Reddening 

Following [PaEerU and lUrrutia et all ()2009f ). we deem 
a quasar "red" if its reddening is E{B — V) > 0.1; 120 of 
the F2M quasars satisfy this criterion. In Figure [15] we 
plot the de-reddened K-hand luminosity for our quasars 
as a function of redshift. We color-code each F2M quasar 
by the amount of reddening, with yellow circles rep- 
resenting lightly reddened quasars, E{B — V) = 0.1, 
and red circles representing heavily-reddened quasars, 
E{B — V) > 1.5. Shades of orange correspond to increas- 
ing amounts of reddening, as annotated in the legend. 
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Fig. 15. — Dereddened X-band absolute magnitude as a func- 
tion of redshift. The colors of the circles correspond to the amount 
of extinction, ranging from low extinction (yellow) to heavily red- 
dened (red). The dotted lines indicate the survey limit {K < 16) 
for increasing amounts of extinction. The small dots are FBQS-II 
and FBQS-III quasars, which we assume are unabsorbed. At every 
redshift, red quasars arc the most luminous. 



We also plot the locations of FBQS quasars (assuming 
their K magnitudes do not experience any extinction) 
with black points. 

Several obvious trends become apparent. At every red- 
shift, F2M red quasars are the most luminous sources. In 
fact, at the most luminous end, the space density of red 
quasars rivals and even exceeds that of blue, unobscured 
quasarS; 

White et al.l ()2003[ ) predicted this result with a sample 
of five red quasars out of 35 FlRST-detected quasars se- 
lected in the /-band. Most of them were normal, unred- 
dened quasars. However, they found five heavily red- 
dened objects at z < 1.3 which, after correcting for red- 
dening, were the most luminous. The interpretation of 
this was that they were "detecting only the most lumi- 
nous tip of the re d quasar iceberg " We began to un - 
cover this trend in iPaper III and in lUrrutia et al.l ()2009l ) 
but these samples were not large enough to trace the full 
range of reddenings across this wide redshift range. We 
discuss the implication of this result in 33 where we in- 
terpret red quasars as a short-lived phase in the context 
of merger-driven QSO ignition and evolution. 

In addition, there appears to be a correlation between 
reddening and redshift in Figure jTSl We see heavily 
reddened objects {E{B — V)> 1.5) at the lowest red- 
shifts, but they become scarce at higher redshifts. This 
is largely a selecti on effect , which has been noted else- 
where (e.g., §6 of IPaper ID . At higher redshifts, what is 
observed as near-infrared light is rest-frame optical emis- 
sion, which is more sensitive to dust extinction. The 
dashed lines in Figure [15] show this explicitly: at a given 
X-band flux limit (e.g., K < 15.5 for 2MASS) in order to 
detect heavily reddened objects at higher redshifts, they 
must be that much more luminous. To reach the red- 
der objects at higher redshifts and overcome this bias, a 
more sensitive near-IR survey is needed. 
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5.3. The Fraction of Quasars That Are Reddened 

We overplot the surface density of F2M quasars as a 
function of their intrinsic K-hand magnitude, after cor- 
recting their apparent K-hand magnitudes, in Figure 
[TUl with open violet squares. We find that after cor- 
recting for extinction F2M quasars comprise 21 ± 2% 
of radio-selected quasars with K < 14.5 magnitudes. 
The 14.5 < K < 15 bin begins to show incompleteness, 
but we can estimate a lower limit to the fraction of red 
quasars with K < 15 at 15±1%. This fraction is re- 
markably consistent with the result from iRichards et al.l 
(gOOS) who find that SDSS misses - 15% of reddened 
quasars that would otherwise be detected in the imaging 
survey. 

This comparison, however, produces a lower limit to 
the fraction of red quasars since the sensitivity of our 
survey to red quasars is strongly luminosity and red- 
shift dependent, as seen by the dashed lines in Figure 
[TSl We therefore should only compare the space density 
of red quasars to FBQS quasars which, when reddened, 
would still be detected in our survey. Our color cuts, 
and the definition we use for "red quasar" , restrict our 
sample to E{B — V) > 0.1. Therefore, we ought not 
count FBQS quasars whose luminosities are too small, 
at a given redshift, to be detected by 2MASS when red- 
dened by E{B — V)=0.1. In other words, we can only 
compare the densities of quasars in the region of AIk 
and z space where the red quasars overlap the FBQS 
points. There may be (and probably are) red quasars at 
lower intrinsic luminosities, but we cannot detect them 
in 2MASS. 

In Figure [16] we plot the surface density as a function 
of Mk of all F2M quasars and compared with FBQS 
quasars that would be detectable in 2MASS if reddened 
by E{B — V)= 0.1. We caution that these curves are not 
luminosity functions since we are plotting a surface den- 
sity (deg~^) not a volume density (Mpc^'^) and we have 
not corrected these number counts by the selection func- 
tions of either survejH. However, since both surveys rely 
on the same flux-limited radio and near-infrared detec- 
tion we assume that they are sufficiently comparable for 
the purposes of measuring the fraction of quasars that are 
reddened. We compare only quasars with Mk < —25.75 
(indicated by the vertical dotted line), since reddened 
quasars with lower intrinsic luminosities would drop out 
of 2MASS. When comparing these two populations, the 
fraction of red quasars is 9.9%. This is, of course, a lower 
limit because of the strong redshift-dependent sensitivity 
of our survey to finding reddened quasars. 

Since our ability to find heavily reddened sources is a 
strong function of redshift, we divide our sample into 
two redshift ranges and compare the distributions of 
quasars with sensitivity limits appropriate to those red- 
shift ranges. In the left panel of Figure [T7] we plot the 
dereddened Mk distributions of red and blue quasars 
with < z < 1.1 and applying the sensitivity limit of 
E{B — V)= 1.0 mag, which Figure [15] shows is valid for 
these lower redshifts. On the right we plot the same dis- 
tributions but for quasars with 1.1 < z < 3.1 and which 
obey the E{B — V)= 0.5 sensitivity limit. Comparing the 
areas under the curves, we find that the surface density of 

We do scale the survey areas by their respective spectroscopic 
completeness. 
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Fig. 16. — Distribution of red quasars as a function of de- 
reddened absolute K-hand magnitude compared with FBQS 
quasars above the 2MASS limit sensitive to E{B — V')=0.1 mag 
(Mfc = —25.75 dotted line). Including all quasars above this limit, 
the fraction of quasars that are reddened is 9.9%. 

the low redshift red quasars is 13.3% of the FBQS coun- 
terpart, while at higher redshifts the surface density of 
red quasars is 19.5% of the FBQS quasars. In both cases, 
and as we noted earlier, we see that red quasars make up 
a higher fraction of quasars at higher luminosities. 

This trend is in direct contradiction to the "receding 
torus" seen in lower luminosity AGN in which a the 
fraction of obscured (and reddened) quasars increases 
with decreasing luminosity. It is hard to envision an 
orientation-based explanation for an increase in the dust 
covering fraction with increasing luminosity. We note 
further that the absence of blue FBQS quasars at the 
very high luminosities seen for F2M quasars is more dif- 
ficult to explain as a selection effect. If the space densities 
of equally luminous unobscured blue quasars is compa- 
rable to or higher than that of red quasars, they would 
be difficult to miss in a blue radio-selected quasar survey. 
We instead interpret the very high de-reddened luminosi- 
ties of F2M quasars in the context of an evolutionary 
picture in Section [7] 

6. BROAD ABSORPTION LINE QUASARS 
Broad absorption line (BAL) quasars show broad 
absorption features (with velocities > 2000 km s^^) 
blueshifted by ~ 3000 - 25000 km s'^ from either C IV 
A1549 and other high ionization UV lines, e.g., N V 
1240A, and/or Mg II A2800 fe.g. IWevmann et al.|[T99l . 
BAL quasars lacking Mg II absorption are classified as 
high-ionization BALs (HiBALs) and those with Mg II 
absorption are low- ionization BALs (LoBALs). An even 
rarer class of LoBALs show absorption fro m metastable 
iron species, Fe II and Fe III (FeLoBALs; [Hazard et al.l 
[19871) . Initially BALs were thought to exist only in radio- 
quiet quasars, although this was likely due to selection 
effects. In radio stacking studies of quasars. [White et al.l 
(j2007) found that the median 20 cm flux density from 
LoBAL q uasars is 2-3 t imes higher than for non-BAL 
quasars. [Becker et al.l (|2000|) found 29 BALs in the 
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Fig. 17. — Distribution of red quasars as a function of de-reddened absolute K-hand magnitude in two redshift bins: (left) < z < 1.1 
above the E(B — V)=l.O sensitivity limit; (right) 1.1 < z < 3.1 above the E{B — V)=0.5 sensitivity limit. The dotted vertical line marks 
the sensitivity limit above the E{B — V) for each plot. We integrate under each curve to this line to compare the space densities of each 
sample. The space density of red quasars in each redshift-reddening bin is 13.3% (left) and 19.5% (right). 



FBQS II survey, amounting to 18% ± 4% of the quasars 
whose redshift range allows their identification. Low- 
ionization BALs are far less abundant than the HiBALs, 
making up between 10% and 20% of BALs depending on 
the selection meth od of their parent q uasar sample (e.g., 
iMenou et alllMlL iTblinp et al.l[2006l ). 

Thirty-three of the F2M red quasars have spectra with 
coverage of Mg II with sufficiently high signal-to-noise 
ratio in the countinuum to measure the presence of broad 
absorption features. This corresponds to objects with 
an optica l spectrum and in th e redshift range 0.9 < z. 
Following lUrrutia et al.l (|2008[ ). we measure the strength 
of the C IV and Mg II ab sorption lines using th e Balnicity 
Index (BI; Appendix A lWevmann et aT]|199lD . a sort of 
modified equivalent width in units of km s~^. 



BI = 



^3000 




/25,000 


0.9 



CdV, 



(4) 



where f{V) is the normalized fiux as a function of veloc- 
ity displacement from the line center, V. This measure- 
ment ignores absorption with velocity widths less than 
2000 km s~^ and may miss so- called mini-BALs . 

To overcome these drawbacks lHall et al.l (|2002l ) defined 
a a more liberal measure of the absorption trough the 
Absorption Index (AI) . We use a modified ver sion of this 
measurement, defined bv lTrump et al.l (|2006[ ): 



AI = - / [1 - fiV)]C'dV, 

'29.000 



(5) 



where f{V) is the normalized template-subtracted flux as 
a function of velocity displacement fro m the line center, 
V (see Appendix A of iHall et al.ll2002l for details on the 
AI). 

In addition, iTrump et al.l (|2006[ ) use a reduced chi- 
squared measurement to determine if an object in their 



sample is a BAL: 



Xtrough 



E- 



1 - fiv) 



(6) 



They require that BALs have Xtrough > 10- 

We compute all three BAL-determining indices for 
both C IV and Mg II for the thirty-three quasars that 
may reveal a BAL. We present these measurements in 
Tabled 

Of these thirty-three quasars, twenty show broad ab- 
sorption lines in their spectra, amounting to 61% of the 
eligible quasars in our sample. They are labeled 'BAL' 
in Column (9) of Table [3l Only one of the BALs found in 
our survey is of the typically more common HiBAL type; 
the rest show absorption in Mg II and are LoBALs. 

Why do we find so many LoBALs in our survey and few 
HiBALs? It has long been known that LoBAL quasars 
have re dder colors than HiBAL qu asars (jWevmann et al.l 
119911 ). ISpravberrv &: Foit3 ()1992[ ) compared composite 
spectra of six LoBALs and 34 HiBALs and found that 
the LoBALs appear to be reddened by an SMC-like dust- 
extinction law with E{B — V) = 0.1. Dust extinction 
preferentially excludes these objects from flux-limited 
quasar surveys and they may be underr epresented by up 
to a factor of ten. 'Becker et al.l ()2000|) find that ~ 3% 
of quasars are LoBALs; ad justing this figure by a factor 
of ten as suggested by iSpravberrv fc Fol"t3 ([1992), im- 
plies that out of thirty-four eligible quasars there may 
be up to 19 LoBALs. This estimate is consistent with 
our result of twenty-one. In fact, the typical extinction 
in LoBALs may be much highe r than that discu s sed in 
ISpravberrv fc Foltd (| 19921) and IWevmann et al.l (|1991D 
since their objects are optically selected. The extinctions 
found for many of these objects from their continuum fits 
suggest much higher reddening than E{B — V) =0.1 and 
may repr esent the large pop ulation of red LoBALs pre- 
dicted bv lBecker et al l 1)2000'). This result, however, does 
not explain the deficit of HiBAL quasars. 
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Since C IV A1549 has a shorter wavelength than Mg II, 
we can only detect HiBALs in z > 1.7 quasars. Of the 
34 quasars able to reveal a BAL, 19 have C IV coverage 
in their spectra. Sixteen of these are LoBALs, one is 
a HiBAL. Many of the remaining objects are so heavily 
reddened that there is no significant flux at or blueward 
of C IV. It is possible that there may be HiBALs among 
these 11 objects, but even so they would still be in the 
minority, far outnumbered by the LoBALs. 

LoBALs and FeLoBALs appear to be fundamen- 
tally different from normal BALs (Beck er et al. 2000; 
White et al.l 120071: lUrrutia et all feoOSt iFarrah et al.1 



20101) TheSEDs of HiBALs do not appear to be different 



from non-BAL quasars suggesting that they are simply 
viewed at an orientation along a disk-wind, which i s ubiq- 
uitou s to all quasars ( White et al. 2000; Gallaghe r et al.l 
I2007D . although they may be slightly reddened by dust 
and gas in the wind. The F2M survey targets red, and 
selects against blue, objects. We may be excluding Hi- 
BALs along with the blue quasars. Since enough contin- 
uum flux must be detected at C IV to detect absorption, 
we both select against HiBALs and have difficulty iden- 
ti fying the ones that a re selected. 

lUrrutia et al.l (|2008D found an anomalously high frac- 
tion (37%) of LoBALs (and FeLoBALs in particular) in 
a sample of FIRST-2MASS-selected red quasars - a sub- 
set of this sample - and suggest these are young objects 
in a "blowout" phase. In our sample, above z ~ 1.7, all 
but one (F2M0134— 0931) of the eligible red quasars are 
BALs, and, with one exception (F2M0854-I-3425), ah are 
LoBALs and most are FeLoBALs. And while some of 
these absorption systems may not be the classical broad 
absorption type (defined by the BI), but they are all dis- 
playing some sort of outflows. In Section [7] we argue 
that this population of objects represents a transitional 
evolutionary stage in quasar evolut ion, which has been 
suggested bv lShen fc Menardl ()20lH ) for quasars exhibit- 
ing associated Mg II absorbers. 



blue, unobscured quasars, also suggesting that they are 
"y oung". Finally, a larg e sample of red quasars reported 
by lUrrutia et al.l ()2008l ) contain an anomalously large 
fraction of LoBALs which they interpret as a "blowout" 
phase representing feedback from the quasar, expelling 
dust and quenching star formation. We find an even 
larger fraction of LoBALs in this, larger, sample of red 
quasars. 

If the red quasars in our survey are a distinct phase 
in quasar evolution then we can place a rough es- 
timate on the duration of the phase in comparison 
to the length of the bl ue quasar phas e , which is es- 
timated^ at ~ 10^ yr ("Y u fc Lul [2005t iMartinil l200l 
Porciani et al.ir2 004: Grazia n et al.l 2004 : iJakobsen et al.l 
20031: iHaimanfc Ceni ,2Q02i: iBaitlik et al.l I1988D . From 



Figure [TU] we integrate the space density of FBQS-II 
Quasars to K = 14.5 magnitudes to be is 6.4 ±0.5 x 10~^ 
deg~^. The extinction-corrected surface density of red 
quasars to K = 14.5 magnitudes is 1.2 ± 0.1 x 10~^ 
deg~^, which is a lower limit since we are not sensitive 
to the most heavily reddened quasars at higher redshifts. 
Nevertheless, the ratio of these space densities suggest 
a rough estimate of the ratio of each phase's duration. 
This is the most direct comparison of quasars in these 
two phases since the FBQS and F2M samples obey the 
same radio and near-infrared fiux limits and differ only in 
the color selection. We find that the ratio of red to blue 
quasars using these samples is 21 ± 2%, implying that 
this emergence phase where the fully obscured quasar is 
s hedding its cocoon last s a few x 10^ yr. 

iHopkins et al.l ()2005b[ ) present a model for the evolu- 
tion of a quasar's light curve during a major merger. 
This model computes the bolometric luminosity, i3-band 
luminosity and column density Nh averaged over many 
lines-of-sight as a function of simulation time. In these 
simulations, during the merger, the quasar bolometric 
luminosity peaks twice. The first peak occurs during a 
heavily obscured phase, where Nh peaks and the B-band 
luminosity is completely attenuated. This completely ob- 



7. RED QUASARS AS AN EVOLUTIONARY PHASE ^^^^^^ phase lasts for ~ 3 - 5 x 10^ years. The second 



Since the discovery of this population of red quasars, 
the evidence suggesting that they are "young" has been 
mounting. Merger-driven models of quasar and galaxy 
co-evolution predict an obscured phase for quasars which 
precedes the famili ar, well-measured, luminous blue 
_guasar phase (e.g., I Sanders et al.l [19881 IHopkins et al.l 
2005b ). T hirteen F 2M red quasars (selected from 



Paper II and IPaper it ) were imaged with HST and re- 



vealed a large fraction (85%) of me rgers (many of whic h 
appear to be major "train wrecks" : lUrrutia et al.ll2008l ). 
consistent with the merger-driven scenario. These thir- 
teen objects have also been observed with Spitzer and 
have low-resolution IRS spectroscopy and MIPS photom- 
etry, which, combined with their spectroscopy, have al- 
lowed a measurement of their black hole masses and bolo- 
metric luminosities and thus estimates of their Eddington 
ratios. These quasars have unusually high Eddington ra- 
tios, and appear to f all below the back hole mass - bulge 
luminosity relation ([Marconi fc H unt 2003.) , indicating 
that they are rapidly growing to "ca tch up" with their 
hosts (Urrutia et al., in preparation). iGeorgakakis et ahl 
(|2009D studied the mid-infrared properties of a small 
sample of red quasars from Spitzer photometry and 
found that they had higher levels of star formation than 



peak in bolometric luminosity corresponds with the max- 
imum observed B-band luminosity and lower Nh, which 
is interpreted as the "normal" , unobscured, blue quasar 
phase. 

The F2M red quasars are not completely obscured, i.e., 
they are not Type 2 quasars; their spectra show broad 
lines (by definition), and some rest- frame i3-band flux 
from the quasar. Their extreme luminosities, disturbed 
morphologies and moderate E{B — V) values suggest that 
they are emerging from the first peak in the quasar's light 
curve. An estimate of the duration of this "emergence" 
phase, based on the Hopki ns et al.l (j2005b[ ) model (e.g.. 
Figure 2), is roughly 5 x 10^ years, consistent with our 
estimate. 

We caution that this merger-driven quasar/galaxy co- 
evolution picture should not necessarily be extended to 
moderate and low luminosity AGN. Recent observations 
of the host galaxies AGN selected by their X-ray emis- 
sion in deep fields whose areas are too small to find 
the rare high luminosity sources presented here (e.g., 
CDFS, GOSMOS) find no evidence for major mergers 
in these (mostly obscured) sources. In fact, the host 
galaxies of these AGN tend to be mostly disk-dominated 
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with some spheroids and the fraction of major mergers 
(~ 1 6%) does not differ from the inactive galaxy p opula- 
tion (iSchawinski et al.ll20lil: iKocevski et al.ll201lD . Fur- 
thermore, low luminosity X-ray selected AGN at z ~ 1, 
which are mostly obscured, have extremely low Edding- 
ton ratios (i/LEdd < 0.1), suggesting that they are in a 
different, "sl ow" phase of growth than the quasars pre- 
sented here (jSimmons et al.|[2011[ ). 

8. SUMMARY AND CONCLUSIONS 

We have identified a complete sample of 120 radio- 
detected, infrared-bright, dust-reddened quasars over the 
9033 deg^ area of the FIRST radio survey. The quasars 
are selected from a candidate hst of 395 FIRST-2MASS 
matches with optical counterparts from the digitized 
POSS-II catalog whose colors obey R — K > 4 and 
J — K > 1.7. Objects that lack an optical counterpart 
are also included and the POSS-II i?-magnitude limit of 
20.80 is used. Reddened quasars are defined as as objects 
with a spectrum that reveals at least one broad emission 
line whose width implies v > 1000 km s""'^ as well as red- 
dening of E{B — V) > 0.1 based on fits to a reddened 
quasar template spectrum. 

The sample spans a broad range of redshifts, 0.13 < 
z < 3.1, and reaches reddenings as high as E{B — V) ^ 
1.55. Compared to radio-plus-optical-selected quasar 
samples such as FBQS and radio-detected quasars in 
SDSS, we find that red quasars make up 11 ± 2% of 
the apparent iiT-band-selected quasars. However, once 
we correct for extinction, we find that, depending on 
how the parent population is defined, red quasars make 
up '~ 15 — 20% of the radio-emitting, luminous quasar 
population. We also find that, at every redshift, red 
quasars are the most intrinsically luminous objects sug- 
gesting that they are in a state of high accretion. This 
will be tested in a future paper where we will estimate 
black hole masses and accretion rates for these quasars. 
We also reproduce the result from Urrutia et al. (2008) 
that the red quasar population contains a large fraction 
of LoBALs, which, evidence suggests, may be young ob- 
jects with strong outflows. 

We therefore interpret dust-reddened quasars as a 
brief evolutionary phase that traces the transition 
from a heavily enshrouded ULIRG-like phase of black 
hole growth to the blue, unobscured quasars found in 
optically-selected samples. If this red quasar population 
is interpreted as an evolutionary phase in the lifetime of 
a quasar, then based on their fraction we estimate the 
red quasar phase to be 15 — 20% of the luminous blue 
quasar lifetime, or a few million years. 
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F2M Red Quasar Candidates 
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(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(1 


00:23:47.96 


-01:45:42.7 


21.44 


19.71 


17.25 


16.05 


15.07 


1.09 


0.79 


2.18 


4.64 


0.342 


Galaxy 


OlO 




00:25:34.91 


-00:52:51.6 


20.90 


19.68 


17.12 


16.56 


15.31 


3.11 


3.47 


1.82 


4.37 




7 


022 


IR 


00:30:09.10 


-00:27:44.4 


19.68 


18.36 


16.50 


15.34 


14.20 


2.95 


2.52 


2.31 


4.16 


0.242 


Starb?ULIRG 






00:31:32.62 


+01:30:01.1 


21.84 


19.79 


17.52 


16.79 


15.75 


1.03 


0.55 


1.77 


4.04 


0.38 


Galaxy 
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15.10 


13.63 
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0.82 


2.94 


5.90 


0.294 


QSO 


05 
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17.47 


16.45 


15.62 


1.11 


1.55 


1.85 


4.27 


0.579 


Galaxy 


022 




00:44:02.81 


-10:54:18.9 


22.34 


19.10 


17.68 


16.07 


15.07 


38.60 


42.57 


2.61 


4.03 


0.431 


NLAGN 
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TABLE 2 

Line Diagnostics of Narrow-Line Spectra 



Object 


Log([0 IIl]/H/3) 


Log([N II] /Ha) 


Log([S II] /Ha) 


Log([0 I] /Ha) 


Diagram 1 


Diagram 2 


Diagram 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


F2M0044-1054 


0.95 




-0.31 


-1.34 




agn 


agn 


F2M0730+2538 


0.73 


-0.68 


-1.02 


-1.91 


agn 


starb 


starb 


F2M0738+2141 


1.31 


-0.29 


-0.61 


-1.18 


agn 


agn 


agn 


F2M0738+3156 


1.04 


-0.41 


-0.86 


-1.53 


agn 


agn 


agn 


F2M0809+6407 


0.68 


-0.50 


-0.87 


-1.38 


agn 


starb 


agn 


F2M0810+1603 


-0.02 


-0.22 


-0.56 


-1.31 


starb 


starb 


starb 


F2M0817+1958 


-1.10 


-1.07 


-1.60 


-1.55 


starb 


starb 


starb 


F2M0818+1305 


1.00 


-0.35 


-0.54 


-1.16 


agn 


agn 


agn 


F2M0820+4728 


-0.29 


-0.13 


-0.45 


-1.12 


starb 


starb 


starb 


F2M0828-I-6420 


0.86 


0.17 


-0.28 


-1.14 


agn 


agn 


agn 


F2M0908+2235'' 


0.77 














F2M0915+1819 


-0.16 


-0.39 


-0.47 


-1.28 


starb 


starb 


starb 


F2M0915+2850 


1.42 


0.23 


-0.08 


-0.83 


agn 


agn 


agn 


F2M0917+2459 


1.08 


-0.27 


-0.31 


-1.08 


agn 


agn 


agn 


F2M0919+1531 


-0.09 


-0.33 


-0.44 


-1.32 


starb 


starb 


starb 


F2M0929-I-6350 


0.16 


-0.63 


-0.39 


-1.50 


starb 


starb 


starb 


F2M0949+2325 


1.76 


—0.24 


-0.31 


-1.33 


agn 


agn 


agn 


F2M1004 +2943 


0.74 


-0.38 


-0.77 


— 1.49 


agn 


agn 


agn 


i< 2M1U15+12U7 


0.98 


—0.46 


—0.28 


—0.88 


agn 


agn 


agn 


F2M1018+2135 


0.50 


— 1.24 


— 1.09 


— 1.58 


starb 


starb 


starb 


F2M1022+1628^ 


0.98 








agn 


agn 


agn 


F2M1049+4010'' 


0.75 














F2M1052-I-0650 


1.99 




-1.17 


-1.05 


agn 


agn 


agn 


F2M1125+0754 


1.12 


-0.65 


-0.50 


-1.42 


agn 


agn 


agn 


F2M1 130+2548 


-0.06 


-1.01 


-1.13 


-1.82 


starb 


starb 


starb 


F2M1139+2743 


0.47 


-0.56 


-0.41 


-0.82 


starb 


agn 


agn 


b 2M1155+08i9 


2.17 




1.56 


2.84 


agn 


agn 


agn 


F2M1201m0332'' 


0.89 














F2M1202+2615'^ 


1.10 








agn 


agn 


agn 


F2M1206p2857 




-1.36 


-1.64 




Starb 


Starb 




F2M1331-0232 


1.14 


-0.27 


-0.49 


-1.03 


agn 


agn 


agn 


F2M1354+0011 


0.55 




-1.33 


-1.80 




starb 


starb 


F2M1446-0050 


-2.54 


-0.23 


0.60 


-0.54 


starb 


starb 


starb 


F2M1448+3056 


0.58 


-0.03 


-0.26 


-0.71 


agn 


agn 


agn 


F2M1518+2336 


-0.01 


-0.41 


-0.62 


-1.71 


starb 


starb 


starb 


F2M1612+3850'' 


1.31 








agn 


agn 


agn 


F2M1638+1058 


0.71 


-0.33 


-0.63 


-1.02 


agn 


agn 


agn 


F2M2328-1107 


0.63 


-0.11 


-0.45 


-0.99 


agn 


agn 


agn 



Note. — Diagram 1 is the classification determined from the [N II] /Ho versus [O IIl]/H/3 diagnostic. 
Diagram 2 is the classification determined from the [S II] /Ha versus [O IIl]/H/3 diagnostic. Diagram 3 is the classification 
determined from the [O I] /Ha versus [O IIl]/H/3 diagnostic. 

^ These objects were at redshifts higher than 0.55 and Ha was outside our spectroscopic range. Nevertheless, we classify these 

sources as AGN since their [O IIl]/H/3 ratios were above the extreme starburst dividing line shown in Figure [7] 

^ These objects are similar to those in except their [O IIl]/H/3 ratios were too ambiguous to determine a clear classification. 



TABLE 3 

Extinction Parameters for F2M Quasars 



Name 




z 


Ks 


Ak 


R-Ks 


E(B - V) 


Spectrum 








(mag) 


(mag) 


(mag) 


(mag) 


(Used) 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


F2M003659.85- 


-011332.3 


0.294 


13.63 


0.40 


5.90 


0.79 


Comb. 


F2M013435.68- 


■093103.0 


2.220 


13.58 


1.30 


>7.22 


0.57 


Comb. 


F2M013613.46- 


-005223.3 


0.718 


15.51 


0.51 


4.64 


0.64 


Comb. 


F2M015647.61- 


-005807.4 


0.507 


14.87 


0.30 


5.01 


0.47 


Opt. 


F2M022950.66- 


-084235.5 


0.440 


14.86 


0.09 


4.32 


0.15 


Opt. 


F2M072910.354 


-333634.0 


0.957 


14.52 


0.83 


5.64 


0.84 


Opt. 
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TABLE 3 

QSO Counts and Surface Densities in K Band 



K Magnitude Range 
(mag) 


FBQS II 
(deg~^ 0.5 mag~i) 


FBQS III 
(deg~^ 0.5 mag~^) 


SDSS^' 
(deg~^ 0.5 mag~^) 


F2M 

(deg~^ 0.5 mag~i) 


Absorption Corrected F2M 
(deg~^ 0.5 mag~^) 


11.0 - 11.5 
11.5 - 12.0 
12.0 - 12.5 
12.5 - 13.0 
13.0-13.5 


0.0012 ± 0.0007 
0.0004 ± 0.0004 
0.007 ±0.001 
0.004 ±0.001 
0.00655 ± 0.002 


0.002 ± 0.002 


0.0003 ± 0.0002 
0.0003 ± 0.0002 
0.0031 ± 0.0007 
0.0053 ± 0.0009 
0.011 ±0.001 


0.0003 ± 0.0002 
0.0001 ± 0.0001 
0.0001 ± 0.0001 
0.0006 ± 0.0003 


0.0004 ± 0.0002 
0.0001 ± 0.0001 
0.0006 ± 0.0003 
0.0023 ± 0.0006 


13.5 - 14.0 


0.0152 ± 0.002 


0.008 ± 0.004 


0.018 ± 0.002 


0.00221°;°°°^ 


0.0033 ± 0.0007 


14.0 - 14.5 


0.033 ± 0.004 


0.036 ± 0.008 


0.038 ± 0.002 


0.0028l°;°°°« 


0.0046 ± 0.0008 


14.5 - 15.0 


0.049 ± 0.004 


0.05 ± 0.01 


0.068 ± 0.003 


0046+° °°'"' 

U.UU4D_Q QQQg 


0.0045 ± 0.0008 


15.0 - 15.5 


0.060 ± 0.005 


0.08 ±0.01 


0.123 ±0.004 


n 007+°-°°^ 


0.0013 ± 0.0004 


15.5 - 16.0 


0.040 ± 0.004 


0.10 ±0.02 


0.124 ±0.005 


0003+° °°^^ 

U.UUUJ_Q QQg2 




-14.5 


0.064 ± 0.005 


0.046 ± 0.009 


0.076 ± 0.003 


0061+° ""^^ 

U.UUDi_Q QQQ9 


0.011 ±0.001 


-15.0 


0.113 ±0.007 


0.10 ±0.01 


0.144 ±0.005 


U.U11_Q QQ^ 


0.016 ± 0.002 



This sample is made up of FIRST-detected quasars from the SDSS. See §4 for how we derive the spax;e densities for these quasars. 



TABLE 4 

Absorption Properties of F2M Quasars with z > 0.9 



Balnicity Index Absorption Index 



Name 


Redshift 


C IV 


Mg II 


C IV 


Mg II 


CIV 


Mg II 


Comment 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


F2M0134-0931 


2.220 


1119 





6112 





0.199 





non-BAL 


F2M0729+3336 


0.957 









390 




4.419 


LoBAL? 


F2M0738+2750 


1.985 








402 


778 


0.025 


83.78 


mini-BAL/LoBAL 


F2M0832+0509 


1.070 

















non-BAL 


F2M0834+5112 


2.391 


798 




2728 




0.249 




FeLoBAL" 


F2M0854+3425 


3.050 


535 




2077 




11.65 




HiBAL 


F2M0904-0145 


1.005 



















F2M0916-0319 


1.560 

















non-BAL 


F2M0916+3854 


1.265 

















non-BAL 


F2M0921+1918 


1.791 


958 





5203 


371 


0.357 


49.59 


miniBAL/FcLoBAL 


F2M0925+4217 


1.879 


232 





2312 


1033 


0.977 


31.36 


miniBAL/FcLoBAL 


F2M0927+3930 


1.170 









1320 




143.0 


miniBAL/FcLoBAL'' 


F2M0943+5417 


2.232 


593 


2539 


4699 


8467 


1.90 


22310 


FeLoBAL 


F2M1004+1229 


2.658 


6121 


11540 


8711 


17374 


4.49 


3062 


FeLoBAL 


F2M1005+6357 


1.280 




6351 




7309 




4.512 




F2M1012+2825 


0.937 




428 




452 




0.611 


non-BAL 


F2M1025-0739 


2.340 








775 





1.369 





non-BAL'^ 


F2M1033+6051 


1.401 

















non-BAL 


F2M1036+2828 


1.762 


541 





2976 


363 


3.921 


1241 


mini-BAL/FeLoBAL 


F2M1222+4223 


2.108 


1409 





2652 


1427 


33.93 


13.48 


mini-B AL /FeLoBAL 


F2M1310-0030 


2.650 


8434 




10752 




19.31 




FeLoBAL 


F2M1311+5513 


0.926 

















non-BAL 


F2M1341+3301 


1.715 





1725 


979 


6368 


0.729 


1031 


FeLoBAL 


F2M1353+3657 


1.311 




4503 




3111 




1.166 




F2M1427+3723 


2.168 


117 





1507 


402 


33.47 


1519 


FeLoBALd 


F2M1456+0114 


2.363 


3991 




6422 




5.84 




FeLoBAL^- 


F2M1500+0231 


1.516 




355 




4232 




10.61 


FeLoBAL 


F2M1531+2423 


2.287 


46 





1383 


351 


13.25 


24.41 


mini-B AL/LoBAL 


F2M1548+0913 


1.509 

















non-BAL 


F2M1549+1245 


2.373 


1604 





4763 


778 


2098 


176.2 


FeLoBAL^ 


F2M1638+2707 


1.689 








49 


579 


0.011 


217.1 


mini-B AL /FeLoBAL 


F2M1650+3242 


1.059 

















non-BAL 


F2M2222-0202 


2.252 








1054 


1184 


69.36 


144.7 


mini-B AL / (Lo)B AL 


■ Classification is 


based on 


Fe II absorption. 













^ Absorption occurs at the quasar redshift, with no velocity shift. 

Three narrow FeLoBAL lines, 
d Mg II mini-BAL; clear C IV BAL, Fe II absorption, 
text 
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76-003007.4 z = 2.650 
:.2 





F2M155102.81 +084401.5 z = 2.507 



E(B-y) = 



0.4 



F2M 1652 14.0 1+443530.6 z = 2.500 



E(B-y) 



0.1 



F2MG834i3.88+51 1214.6 z = 2.39' 
E(B;-V) f 0.4 





0.5 0.6 0.7 0.8 0.9 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

Wavelength (jitm) 



Fig. 6.1. — F2M candidates identified as red quasars, having at least one broad emission line and E{B — V)> 0.1, ordered by redshift. 
The dotted lines show expected positions (in /im) of prominent emission lines in the optical and near-infrared: Lya 1216, N V 1240, 
Si IV 1400, C IV 1550, C III] 1909, Mg II 2800, [O II] 3727, KS 4102, H7 4341, 4862, [O III] 4959, [O III] 5007, Ha 6563, Pa7 10941, 
Pa/3 12822 and Paa 18756A. 
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ii: F2Mil549 

^. E(a-v) 
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~. F2M;1456 
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\[ E(B-V) 
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1 
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F2M153i150.47+242317.6 z = 2.28? 
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Fig. 6.2. — Continued. Spectra of F2M quasars. 
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Fig. 6.3. — Continued. Spectra of F2M quasars. 
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Fig. 6.4. — Continued. Spectra of F2M quasars. 
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Fig. 6.5. — Continued. Spectra of F2M quasars. 
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Fig. 6.6. — Continued. Spectra of F2M quasars. 
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Fig. 6.7. — Continued. Spectra of F2M quasars. 
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Fig. 6.8. — Continued. Spectra of F2M quasars. 
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Fig. 6.9. — Continued. Spectra of F2M quasars. 
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Fig. 6.10. — Continued. Spectra of F2M quasars. 
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Fig. 6.11. — Continued. Spectra of F2M quasars. 
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Fig. 6.12. — Continued. Spectra of F2M quasars. 
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Fig. 6.13. — Continued. Spectra of F2M quasars. 
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Fig. 6.14. — Continued. Spectra of F2M quasars. 
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Fig. 6.15. — Continued. Spectra of F2M quasars. 
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Fig. 6.16. — Continued. Spectra of F2M quasars. 
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Fig. 6.17. — Continued. Spectra of F2M quasars. 
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Fig. 6.18. — Continued. Spectra of F2M quasars. 
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Fig. 6.19. — Continued. Spectra of F2M quasars. 
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Fig. 6.20. — Continued. Spectra of F2M quasars. 



